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Progress in the theory of anomalous diffusion in weakly turbulent cold magnetized plasmas is 
explained. Several proposed models advanced in the literature are discussed. Emphasis is put on 
a new proposed mechanism for anomalous diffusion transport mechanism based on the coupled 
■ action of conductive walls (excluding electrodes) bounding the plasma drain current (edge diffusion) 

together with the magnetic field flux "cutting" the area traced by the charged particles in their 
orbital motion. The same reasoning is shown to apply to the plasma core anomalous diffusion. 
£NJ ■ The proposed mechanism is expected to be valid in regimes when plasma diffusion scales as Bohm 

£H ' diffusion and at high B/N, when collisions are of secondary importance. 
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r-j 1 I. INTRODUCTION 

Any phenomena occurring with interfacial systems has a fundamental importance in science and technology [l[ (e.g. 
electrostatic charging of insulators, surface tension, forward conduction in p-n junctions). Specifically, the problem of 
' the plasma-wall interactions is of major importance in plasma physics. 

Historically, the anomalously high diffusion of ions across magnetic field lines in Calutron ion sources (electro- 
magnetic separator used by E. Lawrence for uranium isotopes) gave the firsts indications of the onset of a new 



mechanism [2| . It has been noticed that the plasma moves across the magnetic confining field at a much higher aver- 
. i-h ' age velocity than it is predicted by classical considerations. The classical diffusion coefficient is given by D± = r/p/B 2 , 
while the anomalous diffusion coefficient by D± = akT/B. Manifestly it is needed a better understanding of the 
^ (— i physical laws governing matter in the far nonequilibrium state, and this is a challenging issue for the advancement of 
pi ' this frontier of physics. 

'— 1 . General proposed explanations were advanced. The first one was Simon's " short-circuit" problem, suggesting that 
the observed losses could be explained by the highly anisotropic medium induced by the magnetic field lines, favoring 
electron current to the conducting walls 0. Experiments done by Geissler Q in the 1960's have shown that diffusion 
1 in a plasma across a magnetic field was nearly classical (standard) diffusion when insulating walls impose plasma 
\ ambipolarity, but in the presence of conducting walls charged particles diffused at a much higher rate. 

This problem of plasma-wall interaction becomes more complex when a complete description is aimed of a magne- 
tized nonisothermal plasma transport in a conducting vessel. Beilinson et al. [5[ have shown the possibility to control 
, the discharge parameters by applying potential difference to sectioned vessel conducting walls. 

In the area of fusion reactors, there is strong indication that for plasmas large but finite Bohm- like diffusion 
ON . coefficient appears above a certain range of B Q . Experiments give evidence of transport of particles and energy to 
the walls]?]] • At the end of the 1960s, experimental results obtained in weakly ionized plasma [1] and in a hot electron 
plasma [8j (this one proposing a possible mechanism of flute instability) indicated a strong influence that conducting 
walls have on plasma losses across magnetic field lines. Geissler [J] suggested that the most probable explanation was 
the existence of diffusion-driven current flow through the plasma to the walls. Concerning fusion reactions, Taylor Q 
provided a new interpretation of tokamak fluctuations as due to an inward particle flux resulting from the onset of 
filamentary currents. 

Progress in the understanding of the generation of confinement states in a plasma is fundamental [10( to pursue the 
dream of a fusion reactor [TTI.Il2||. Anomalous diffusion is a cornerstone in this quest, as recent research with tokamaks 
suggest that the containment time is r ~ 10 8 i? 2 /2Db, with R denoting the minor radius of a tokamak plasma and 
Db is the Bohm diffusion coefficient 13]. Controlled nuclear fusion experiments have shown that transport of energy 
and particles across magnetic field lines is anomalously large (i.e, not predicted by classical collision theory). 
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The conjecture made by Bohm is that the diffusion coefficient is Db = akT/eB, where T is the plasma temperature 
and a is a numerical coefficient, empirically taken to be 1/16 Q. Initially, the origin of the anomalous diffusion was 
assumed to be due to the turbulence of small-scale instabilities (see, for example, Refs. [fiHH, EH). However, it is now 
clear that there is a number of different mechanisms that can lead to anomalous diffusion such as, coherent structures, 
avalanches type processes and streamers, which have a different character than a purely diffusive transport process. 
Recent experimental results such as scaling of the confinement time in L-mode plasmas and perturbative experiments 
undermine the previous paradigm built on the standard transport processes [la . [l6j showing conclusively that there 
are many regimes where plasma diffusion does not scale as B~ x . 

This paper put emphasis on a mechanism of wall current drain set up together with the magnetic field "cutting" 
lines across the area traced by the charged particles trajectories. The proposed mechanism of anomalous diffusion 
is expected to be valid in purely diffusive regimes when plasma diffusion scales as Bohm diffusion, both in the edge 
and core of a cold magnetized plasma. At his stage it was considered of secondary importance the role of collisions 
in randomizing the particle's distribution function. From collisional low temperature plasmas to a burning fusion 
plasma subject the plasma confinement vessel to strong wall load, both in stellarator or tokamak operating modes, 
this explanation could be of considerable interest, particularly when diffusive transport process are dominant. 

We would like to stress that this work is not free from omission of important contributions. 



II. SIMON'S "SHORT-CIRCUIT" THEORY 



The first attempt to explain why the plasma diffuse at a much higher average velocity than it is predicted by 
classical theory has been advanced by A. Simon j3(. The magnetic field lines structure a highly anisotropic medium. 
Any fluctuation of the space charge builds up an electric field, which has a strong effect on the currents parallel to 
the magnetic field lines. In fact, the classical equation for conductivity across a magnetic field is given by 

1 + -k 

where Q c = eB/m is the electron cyclotron frequency, v e is the electron collision frequency, and <tq = e 2 n e /mi' e is 
the conductivity in the absence of a magnetic field. By the contrary, due to Sl c /v e S> 1, this electric field is too 
small to have any importance on the crossed-field conductivity. From this results that there is a strong current to 
the wall without a concomitant current to different regions of the plasma, making of this situation a kind of circuital 
"short-circuit" problem. Although Simon attempted to explain the anomalously high rate of diffusion in Calutron 
ion sources in the frame of the classical diffusion theory calculating the coefficient D± as being approximately equal 
to the transverse diffusion coefficient of the ions. His proposal is not suitable, however, because the experimental 
determination of D± by means of a decaying plasma have shown that according to the magnetic field strengthen 
the transverse diffusion coefficient can be much higher than the classical one or smaller than the transverse diffusion 
coefficient of the ions H . 



III. PLASMA TURBULENCE AND TRANSPORT 



Purely diffusive transport models cannot give convincing explanations for a variety of experiments in magnetically 
confined plasmas in fusion engineering devices, particularly the scaling of the confinement time in L-mode plasmas. 
The assumed underlying instabilities are driven by either the pressure gradient or the ion temperature gradient. It 
is well established fact that transport in high temperature confined plasmas is driven by turbulence and plasma 
profiles, and are subject to transition from L-mode to H-mode (characterized by a a very steep gradient near the 
plasma surface) fl7| . The non-linearity in the gradient-flux relation is the source of turbulence and turbulence-driven 
transport. The fluxes contain all the dynamic information on the transport process. Accordingly, changes in the 
gradient trigger local instabilities in the plasma. This local instability induces an increase in the nearby gradients, 
thus causing a propagation of the instability all across the plasma. In particular, an excessive pressure on the core 
propagates to the edge in a kind of avalanche. 

In weakly turbulent cold magnetized plasmas, besides the Calutron ion sources and the magnetron, the study of 
particle transport in crossed electric and magnetic fields results from applications to electromagnetic space propulsion 
(Hall thrusters). Those plasma accelerators work with a radial magnetic field that prevents electron flow toward the 
anode and forcing the electrons closed-loop drift around the axis of the annular geometry. Neutrals coming from 
the anode are ionized in this rotating electrons cloud, while ions are accelerated by an axial electric field that freely 
accelerate them out from the device. This effect develops in the so called extended acceleration zone (or electric- 
magnetic region plasma). In this acceleration zone the electron gyro-radius and the Debye shielding length are small 



3 



relative to the apparatus dimension, while the ion gyro radius is larger than the apparatus typical length. From these 
spatial scales results that the electron motions are [E x B] drifts, but ions are accelerated by the electric field that 
develops in the plasma. The first observations of a large amount of electron transport toward the anode have been 
noticed in the 60's (e.g., Ref. [HI]) and they have been related to electric field fluctuations since they were correlated 
with the density variations in order to produce anomalous transport. Other possible mechanism that could explain 
the high electrons transverse conductivity were advanced: collisions with the wall nam. Electrons moving freely 
along the lines of forces of the magnetic field collide with the wall more frequently than with ions and neutrals, being 
reflected at the wall and enhancing emission of low-energy secondary electrons from the wall. As referred, the other 
strong candidate, which could possibly be the source of a higher axial electron current than predicted by the standard 
classical kinetic theory, is the turbulent plasma fluctuations. But it seems that there is no clear consensus on this 

issue mmm. 

The magnetron is a sputtering tool, used for reactive deposition and etching. The magnetron effect is applicable 
to different geometries and only need a closed-loop [E x B] drift to work. Rossnagel et al. [24[ have shown that the 
Hall-to-discharge current ratio measured in those configurations could be explained if the high collision frequencies 
for electrons were associated to Bohm diffusion. In particular, Kaufman [25| argues that anomalous diffusion in 
closed-loop [E x B] thrusters could shift from core diffusion to edge diffusion (or wall effects) with increasing magnetic 
fields. 



IV. CIRCUITAL MODEL OF ANOMALOUS DIFFUSION 

In a seminal paper [26| a conjecture was proposed based on the principle of minimum entropy-production rate, 
stating that a plasma will be more stable whenever the internal product of the current density j by an elementary 
conducting area dA at every point of the boundary - excluding the surface collecting the driving current - is null, 
(j • A) = at any point of the boundary (and excluding the surfaces collecting the discharge current), independently 
of the resistance Ri. The general idea proposed by Robertson (26[ assumes that the plasma boundary is composed of 
small elements of area Ai, each one isolated from the others, but each one connected to the exterior common circuit 
through its own resistor Ri and voltage Vi. The entropy production rate in the external circuits is given by: 

f =E?(ji-^) 2 ^, (2) 

i 

where T is the temperature of the resistors, supposed to be in thermal equilibrium with all the others. It is important 
to remark that the summation is over the different conducting areas eventually confining the plasma, excluding the 
electrodes areas. FigQ] illustrates this concept. 

We consider a simple axysymmetric magnetic configuration with magnetic field lines parallel to z-axis with a plasma 
confined between two electrodes (see Fig.l). In general terms, a particle's motion in the plasma results in a massive 
flux. As long as the flux is installed, the flux will depend naturally on a force F - in this case the pressure gradient- 
driven process of diffusion to the wall - responsible of the wall driven current j. According to the fundamental 
thermodynamic relation, the plasma internal energy variation dU is related to the amount of entropy supplied or 
rejected and the work done by the driven force, through the equation: 

d 4 = li .AfR + ( F - d 4)- (3) 



dt V dt. 

The last term we identify with the macroscopic diffusion velocity depicting the process of plasma expansion to 
the wall. To simplify somehow the calculations we assume a single plasma fluid under the action of a pressure gradient 
(F = ALdp/dy j , where j is a unit vector directed along the Oy axis). 

In the presence of steady and uniform magnetic field lines (this simplifies the equations, but do not limit the 
applicability of the model), the particles stream freely along them. From magnetohydrodynamics we have a kind of 
generalized Ohm's law (see, for example, Ref. (27j): 

Vp = -enE - en[v x B] + [j x B] - — , (4) 

a 

where a = e 2 nr e /m e is the electric conductivity, with r e denoting the average collision time between electrons and 
ions. The force balance equation is given by: 



Vp= [j x B], 



(5) 
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FIG. 1: Schematic diagram of the plasma boundary connected to the common circuit through conducting walls. 



valid whenever the Larmor radius is smaller than the Debye radius. This assumption simplifies further the extension 
of our model to high enough magnetic fields. Therefore, after inserting Eq. [5] into Eq. Q]the y component of velocity 
is obtained: 

E x 1 dp 

V ^-B'^d^- (6) 

From Eq.[5]we have the classical diffusion coefficient scaling with 1/B 2 and thus implying a random walk of step length 
Tl (Larmor radius). To get the anomalous diffusion coefficient and as well understand better its related physics, we 
must consider the process of diffusion to the wall - in the presence of an entropy source - with the combined action of 
the wall current drain, as already introduced in Eq. \3\ 

Therefore, using the guiding center plasma model the particle motion is made with velocity given by: 

[Vg x B] 

J = env d = . (7) 

This equation forms the base of a simplified theory of magnetic confinement. In fact, the validity of Eq.[7]is restrained 
to the high magnetic field limit, when the Larmor radius is shorter than the Debye radius. 

Considering the motion along only one direction perpendicular to the wall (y-axis), it is clear that 

If we consider a quasi-steady state plasma operation, the plasma total energy should be sustained. Hence, dU / dt = 0, 
and the power associated with the driven pressure-gradient is just maintaining the dissipative process of plasma losses 
on the wall. Eq. [3] governs the evolution of the diffusion velocity. Hence, we have 

nRA kT dn dn 

with Dt denoting the transverse (across the magnetic field) diffusion coefficient given by: 

nRA kT , s 

D T = — w (10) 
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This new result coincides with the classical diffusion coefficient [28[ whenever nRA/L = rav^je 2 , containing a 
dependence on collision frequency and particle number density. Other theoretical approaches to this problem were 
advanced by Bohm Q , who proposed an empirically-driven diffusion coefficient associating plasma oscillations as the 
source of the enhanced diffusion, while Tonks [2^ | have shown that the current density that is present in a magnetically 
immobilized plasma is only generated by the particle density gradient, not being associated with any drift of matter. 
Simon electron "short-circuit" Q scheme attempt to explain the different rates of diffusion, electrons and ions do 
experiment across the magnetic field. While the ion flux dominates the radial diffusion, the electron flux dominates 
axial losses, due to an unbalance of currents flowing to the wall. 

In the absence of collisions, the guiding centers of charged particles behave as permanently attached to the same 
lines of force. On the contrary, as a result of collisions with others charged particles the guiding centers shift from one 
line of force to another resulting in a diffusion of plasma across the field lines. In our model, each orbit constitutes 
an elementary current I eventually crossing the wall. 

However, the particle diffusion coefficient as shown in Eq. [10] gives evidence of an interplay with the resistance 
that the elementary circuit offer when in contact with the walls in the presence of the frozen-in effect. In fact, for 
sufficiently strong magnetic fields apparently a hydrodynamic behavior of the plasma is installed [H, [33] , with the 
appearance of " convective cells" and the 1/B behavior dominates, giving birth to the anomalous diffusion mechanism. 
The onset of freezing magnetic lines is valid whenever the Lundquist number S > 1 (convection of the magnetic field 
dominated medium). In this case the magnetic field lines are frozen-in in the medium (consequence of a vortex type 
of character of the magnetic field B) and the flux of them across a given surface is constant: 

$ = BA' = BL 2 a. (11) 

Remark that A' is now the surface delimited by the elementary circuit 7 (see Fig. [2]) and a < 1 is just a geometrical 
factor (e.g. a = 7r/4 at the limit of a circular orbit). This situation is fundamental to the onset of anomalous diffusion. 
Free electrons orbits are helical, but as Fig. [2] shows, their projections at rigth angles to the field are circular. Each 
particle orbit constitute an elementary circuit with B-field cutting its surface being associated with it an elementary 
flux $. At the same time we can envisage each orbit as constituting by itself an elementary circuit, some of them 
intersecting the wall and thus the circuit is closed inside the wall. Therefore a resistance R drags the charged flow 
at the conducting wall. It is therefore plausible to associate to this elementary circuit a potential drop V and all the 
process being equivalent to a current I flowing through the elementary circuit. 

Assuming the plasma is a typical weakly coupled, hot diffuse plasma with a plasma parameter (number of particles 
in Debye sphere) A = nX 3 De w 1, it is more likely to expect nearly equal average kinetic and potential energy. However, 
the typical plasma parameter encountered in glow discharges or in nuclear fusion is A ^ 1. This means that the 
average kinetic energy is larger than the average potential energy. To contemplate all range of A we can relate them 
through the relationship 

pV = (J • A)*. (12) 

Here, p is the charge density, A is the vector potential, J is the current density and 5 < 1 is just a parameter 
representing the ratio of potential to kinetic energy. Of course, when A > 1, then 5 < 1. This basic assumption is 
consistent with the hydrodynamic approximation taken in the development of equations. The limitations of the model 
are related with the unknowns A and 5 that can be uncovered only through a self-consistent model of the plasma. 
However, our analysis of anomalous diffusion remains general and added new insight to the phenomena. 

Now suppose that the diffusion current is along y-axis J = — J y u y (see Fig.l). Consequently, A = —A y u y , and 
then the potential drop will depend on x-coordinate: 

p[V{xx) - V(x )} = J y [A y ( Xl ) - A y (x )]S. (13) 

Multiplying both members by the area A' = x\Z\ and length L = y±, we have 

QAV = Iyi[A v { Xl ) - A y (x )}5 = 1*5. (14) 

$ = <f> (A • g?x) is the flux of the magnetic field through the closed surface bounded by the line element chc (elementary 
circuit 7, see also Figf2]). By other side, naturally, the total charge present on the volume V = x\y\Z\ is such as 
Q = ie, with i an integer. This integer must be related to ions charge number. From Eq. [14] we obtain 

AV $ BL 2 , s 

R= — = S- = aS . 15) 

1 Q ie 

But, the particle density is given by n — N/LA, with N being now the total number of charged particles present in 
volume V = AL. Since i — N, we retrieve finally the so-called Bohm-diffusion coefficient 

kT 

D B =a6—. (16) 
eB 
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FIG. 2: Schematic of the geometry for the plasma-wall current drain model. The uniform magnetic field points downward 
along Oz. Particles describe orbits in the plane xOy intersecting the wall (plan xOz). Orbits are represented by a semi-circular 
line for convenience. L is the maximum distance the trajectory attains from the wall. 



So far, our arguments were applied to edge anomalous diffusion. But they can be generalized to the core anomalous 
diffusion processes, provided that diffusive transport processes are dominant. For this purpose consider instead of a 
conducting surface a virtual surface delimiting a given volume, as shown in Fig. [3J 

Our coefficient is time-dependent and can be written under the form: 

kT(t) 9(t) 

d± = s ^bW^- (17) 

The nonrelativistic solutions of dynamical equation of a charged particle in time-dependent but homogeneous electric 
and magnetic field give the following approximative expressions for the width trajectories along the xOy plane (see, 
e.g. Kef. EH) by: 

Ax=^cos(n c t + X ) A2/ = Tf£sin(Q c t + x), (18) 

where x is the initial phase, v± denotes the component of the velocity perpendicular to the magnetic field and the =F 
sign applies to electrons (-) or positive ions (+). From them we can retrieve the flux "cutting" area: 

1 / ft 



AnAx.Ay=-- sin(2fi c t + 2 X ). (19) 

Then the anomalous diffusion coefficient is just given by: 

As we can see in FigUthis last expression describes fairly well the diffusion process for high enough B/N values 
(the magnetic field to gas number density ratio) and explains the main processes building-up such effects as: i) the 



FIG. 3: Volume control and particle's trajectory submitted to magnetic field. Magnetic field lines point downward. 



negative diffusion, which results from the contraction of the flux "cutting" area; ii) the ciclotronic modulation imprint 
on the transverse diffusion coefficient; iii) and the anomalous diffusion, due to the fast flux rate of the magnetic field 
through the area A) . All this signs can be seen on Fig. [J] were it is shown a comparison of numerical results (5000 
Hx, 1 Hx=10~ 27 T.m 3 ) obtained with Monte Carlo simulations of electron transport in crossed magnetic and electric 
fields by Petrovic et al. [32l . l33j with the theoretical prediction given by our Eq. [501 As long as only a self-consistent 
model could give us an exact value of the ratio of potential to kinetic energy <5, we assume here S = 1/40. The 
full agreement with the numerical calculations is not obtained due to neglecting effects related to the electric field 
variation in time and of the assumed collisionless approximation. This explains the big discrepancy shown in Fig. [4] 
when compared with the diffusion coefficient at 1000 Hx when collisions begin to be far more important to randomize 
individual trajectories and our approach is no more valid (at low enough B/N values). 



V. DISCUSSION AND SUMMARY 



However, Eq. 1161 suffers from the indetermination of the geometrical factor a. This factor is related to the ions 
charge number, it depends on the magnetic field magnitude and as well as on the external operating conditions (due 
to increased collisional processes, for ex.). The exact value of the product aS can only be determined through a 
self-consistent plasma model, but we should expect from the above discussion that ad < 1. For a 100-eV plasma in 
a 1-T field, we obtain Db = 1-67 m 2 /s (using the thermal to magnetic energy ratio with particle's density n = 10 14 
cm -3 ). Furthermore, Eq.[l5]can be used as a boundary condition (simulating an electrically floating surface) imposed 
when solving Poisson equation. 

Also it worth to emphasize that when inserting Eq.[T5]into Eq. 1101 and considering the usual definition of momentum 
transfer cross section, then it can be obtained a new expression for the classical diffusion coefficient as a function of 
the ratio of collisional v and cyclotron frequency fi, although (and in contrast with the standard expression), now 
also dependent on the geometrical factor a and energy ratio 5: 

v kT 

D T = (aS)-—. (21) 
il m 

This explains the strong dependence of the classical diffusion coefficient on v/Q showing signs of anomalous diffusion 
as discussed in Ref. [32J (obtained with a time resolved Monte Carlo simulation in an infinite gas under uniform 
fields) and, in addition, the strong oscillations shown up in the calculations of the time dependence of the transverse 
component of the diffusion tensor for electrons in low-temperature rf argon plasma. Those basic features result on one 
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FIG. 4: Comparison between numerical results for 5000 Hx obtained by Monte Carlo simulations Refs. [32,33] of electron 
transport in crossed electric and magnetic rf fields in argon and Eq. [20]-theo. Dashed line: external time-dependent magnetic 
field. Parameters used: 5 — x\ applied frequency, f=100 MHz; ciclotron frequency, Q c ~ 10 10 Hz; ^ = 5.4 eV; p=l Torr, 
T g = 300 K, x=0. 



side from its dependence on R, which is proportional to the flux. Therefore, a flux variation can give an equivalent 
effect to the previously proposed mechanism: whenever a decrease (or increase) in the flux is onset through time 
dependence of electric and magnetic fields, it occurs a strong increase (or decrease) of the diffusion coefficient. By 
other side, when the resistance increases it occurs a related decrease of charged particles tangential velocity and its 
mean energy. So far, this model gives a new insight into the results referred in |32j and also explains why the same 
effect is not obtained from the solution of the non-conservative Boltzmann equation as applied to an oxygen magnetron 
discharge with constant electric and magnetic fields [3^ |. 

A further application of Eq. [5] to a cold plasma can give a new insight into the " ambipolar-like" diffusion processes. 
Considering just one conducting surface (besides the electrodes driving the main current into the plasma) and the 
plasma build-up of electrons and one ion component to simplify matters, we obtain: 

dS e 2 

— = — (-n e /i e E + D e Vn e + mfMiE - D l Vn l ) 2 A 2 R. (22) 
at 1 

Under the usual assumptions of quasi- neutrality and quasi-stationary plasma (see, for example, Ref. [28|), the 
following conditions must be verified: 



(23) 



and hence, Eq. [2] becomes: 



^ = ^[E(eMi - MeK + Vn e (D e - D % e)] 2 A 2 R. (24) 
at 1 

For a stable steady-state plasma with no entropy sources the condition S = prevails and then an " ambipolar-like" 
electric field is recovered [281 ] : 

E= D e -eD l Vn 1 

It means that the conducting surface must be at its floating potential. Such conceptual formulation provides new 
insight into "ambipolar-like" diffusion processes. In a thermal equilibrium state, a plasma confined by insulating 
walls will have an effective coefficient given by the above Eq. 1251 a situation frequently encountered in industrial 
applications. This example by itself relates ambipolar diffusion with no entropy production in the plasma. However, 
allowing plasma currents to the walls, entropy production is greatly enhanced generating altogether instabilities and 
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plasma losses [26j]. As long as confined plasmas are in a far-nonequilibrium state (with external surroundings) it 
is necessary to establish a generalized principle that rule matter, and this circuital model for anomalous diffusion 
represents some progress in the physics of plasmas as nonequilibrium systems. 

To summarize, we introduced in this study a simple circuital mechanism providing an interpretation of the anoma- 
lous diffusion in a magnetized confined plasma in a purely diffusive transport regime. The coupled action of the 
magnetic field "cutting" flux through the areas traced by the charge carriers elementary orbits, together with the 
elementary electric circuit constituted by the charged particle trajectory itself are at the basis of the anomalous dif- 
fusion process. Whenever conducting walls bounding the plasma drain the current (edge diffusion) or, at the plasma 
core, the magnetic field flux through the areas traced by the charged particles varies, a Bohm-like behavior of the 
transverse diffusion coefficient can be expected. Eq. [5D]can be used as an analytical formula when simulating plasma 
behavior at high B/N. In the near future we hope to generalize this model taking into account random collisions. 
The suggested mechanism could lead to a better understanding of the mechanism of plasma-wall interaction and help 
to develop a full-scale numerical modeling of present fusion devices or collisional low-temperature plasmas. 

The author gratefully acknowledge the data supplied to us by Zoran Lj. Petrovic and Zoran Raspopovic used in 
our Fig. 4 . 
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